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Abstract 

We compute the J/ip + 'y invariant-mass distribution from the QCD subprocess 
g + g ^ J /'>p + 7- At large masses, this subprocess is the dominant mechanism 
for J/ijj + ^ production, and data could provide a good test of QCD. The mass 
distribution peaks at relatively small masses (3.4-4.0 GeV) and the subprocess 
could, therefore, represent a significant QCD background to J/'i/'+7 decay of heavier 
charmonia. We also analyze the J/^ angular distribution in the J + ^ rest frame. 
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There has been much interest recently in associated production of a. J/tp and a photon 
(7) in colhder 0] and fixed-target experiments 0, 0]. In the framework of the colour- 
singlet model 1^ 1^, the QCD subprocess that contributes to the production of a J/ip and 
a photon is: 

g + g^J/^ + -f. (1) 

The above subprocess is related to that for photoproduction of J/ip by crossing. The 
photoproduction process and the process in Eq. (|l]) are gluon-initiated, and have been 
suggested as sensitive probes of gluon densities in lepton-nucleon experiments @], and in 
hadron-nucleon experiments [Q, 0, respectively. 

The J/tp + j mode is also the decay mode in which the P-state charmonia are detected 
1^. In this letter, we emphasise that, at low values of the invariant mass of the J/ip + 'y 
pair, the process in Eq. (|l|) forms an important (and computable) part of the background 
to the decay of P-state charmonia into J/ip and 7. We present in the following the 
invariant-mass distributions for the process in Eq. (|l]). In the mass region of the x's, this 
process is a background to the x decay signal. At larger masses the J/ip + •y continuum 
is dominated by the QCD subprocess and could serve as an interesting test of QCD. 
The differential cross-section for the process g + g^J/tp + j is 



dMinvdxrdxi 



(2) 



where Mj„^ is the invariant mass of the J/ip-j pair, xt = 2pt/\/s (a/s being the centre- 
of-mass energy), and G{xi^2) is the gluon distribution in the beam and the target, respec- 
tively. The subprocess cross-section is [H 
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The modulus squared, |P(0)p, of the wave-function at the origin is related to the leptonic 
decay width by 



|P(0)p = 

and J-' in Eq. (Q), is given by 
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In the above equation, yi is the rapidity of the J/ip, r = Mj/s, and xt = yXj- + 4r. The 
rapidity of the J/ip is expressed in terms of the variables Xi, Minv and xt by the following 
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kinematic relation 



yi = In 

with 



[X1X2 + t) + \l[xiX2 + r)2 - XiX2x']^ 
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Using the cross-section in Eq. we obtain the invariant-mass distributions by inte- 
grating over Xi and xy, with the ranges of integration chosen to be 



^ <xi <1; 



<xr <^-^^l^^ 
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Note that the xt integration is not bounded from below, and we are integrating over 
regions of phase space where the photon could become soft. The subprocess cross-section 
in Eq. @ is finite in the soft-photon limit. Nevertheless, for perturbation theory to provide 
a reliable answer, internal lines in the Feynman graphs of @] must remain sufficiently off- 
shell. This requirement means that the kinematic region should be chosen such that the 
photon does not become soft. This condition can be ensured if we choose Mj„„ such that 

ML = Mj + A, (9) 

where A is larger than some value, e.g. 1 GeV^. If M^nv > 3.25 GeV, then the cross- 
section given in Eq. can be considered a reasonable lowest-order estimate of J/'?/' + 7 
production. Contributions from soft-gluon radiation [Q will tend to increase the cross- 
section, and so our estimate can be regarded as a lower bound on the cross section arising 
from the gluon-gluon fusion subprocess. A more reliable estimate of the cross-section 
could be obtained after the soft-gluon corrections to this process are computed. We plan 
to undertake a study of these corrections in the future. 

Integrating x\ and xt over the ranges specified above, we obtain the invariant- mass 
distributions for Mj„^ > 3.15 GeV. We present results for 'pN and vrA^ collisions, for typ- 
ical values of \fs. For the proton structure functions, we use the MRS-SO distributions 



and for the pion distributions, the GRV-P distributions |T0|. The parton distributions 
are evolved to the scale = Mj +p%. In Fig. 1 we show the invariant-mass distribution 
Mf^^da I dMinv as a function of {Minv — Mj)/ y/s, for pN collisions. The two curves corre- 
spond to beam energies of 300 GeV and 600 GeV, respectively. We see that approximate 
scaling sets in, in the region beyond {Mmv — Mj)/y/s ~ 0.07. 

In Fig. 2 we plot the number of events/GeV for pN collisions and compare with data 
from an ISR experiment |^ and the E705 Fermilab experiment 0, {^/s = 62 GeV and 
23.72 GeV, respectively). For the ISR data we display the distribution as a function of 
Minvi whereas for the E705 data we use AM = Mj„^, — Mj. In comparing the shapes of 
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our curves with the data, we normahse to the experimental points at Minv = 4.8 GeV and 
AM = 0.88 GeV, respectively. The shape of the distributions away from the resonance 
region may be in reasonable agreement with the data, but it is desirable to have data 
that extend to larger values of Mj^^. This experimental information will be valuable 
since it will serve as a test of the QCD predictions of both the magnitude and the Mj„„ 
dependence. We suggest that data on associated production of J/V' + 7 at large invariant 
masses could provide a QCD test of gluon initiated hard-scattering reactions similar in 
quality to the QCD test that massive lepton-pair production provides for quark-antiquark 
initiated hard-scattering. 

With the normalisation fixed as stated in the paragraph above, we find that the QCD 
process could account for about 25% of the events under the resonance, for both sets of 
data. In 1^ the background to x decay is constructed from uncorrected J/ip and 7 
momentum vectors obtained from different events and normalized to the data above the 
X mass region. Thus, in the analyses of |^, the QCD signal that we are discussing 
seems to have been subsumed into the "background". 

In Fig. 3 we show the invariant- mass distributions for both pN and vrA^ collisions at 
y/s = 23.72 GeV. A lower cut on the photon energy is often employed in the data because 
definition of the photon may be poor at low photon energies. We present results for three 
different cuts on the laboratory energy of the photon {e}^^^ > 0, 1.0, 2.5 GeV). The mass 
distribution is not very sensitive to the photon energy cut, although the peak is shifted 
slightly to larger Mj„^ as the photon energy cut is raised. 

In addition to the mass distribution discussed above, we have analysed the angular 
distributions of the J/ip in the J/ip-^ rest-frame (which corresponds, in lowest order 
perturbation theory, to the subprocess centre-of-mass frame) . The angular distribution is 
expressed as 

_ (5-^^3)^^G'(xi,g2)rf(x2,Q2)^, (10) 



dcos9*dMinvdxi XiMinv di 
with pt given in terms of the polar angle 6* of the J/ as 

,,= (£^,i„9.. (u) 

2M- 

After integrating over xi, we obtain da/dcos6*dMinv In Fig. 4 we show the angular 
distribution as a function of cos6'* for different values of Mj„^, for pp coUisions at i/i = 
23.72 GeV. The distribution is symmetric about cosO* = 0, and the dependence on cos^* 
becomes stonger with increasing Mj„^. 

Experimental information on angular distributions in the J/ip-j rest system should 
be useful in establishing the production mechanism. In the decay of X2's to J/ip + 7, the 
J/ip angular distribution is expected to follow ||1T| a (1 + cos^6'*) behaviour. In the same 
mass region (Mj^^, = 3.5 GeV), the QCD angular distribution that we have computed 
obeys a (1 + acos^^*) fit, with a = 0.48. Selections on cos6'* should help to enhance the 
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resonance signal relative to the QCD process and vice-versa. For larger Mj„„, the QCD 
angular distribution shows a large contribution from an additional cos^^* term. In the 
X resonance region, the resonance and part of the QCD background occur in the same 
partial wave. A complete analysis of this mass region would require a proper treatment 
of the effects of final-state interactions. [^] 

In summary, we have presented calculations of invariant- mass distributions for asso- 
ciated production of J/^ + 7 through gluon-gluon fusion. This QCD process provides a 
potentially important background to the J/V' + 7 decay mode of the P-state charmonia. 
Angular distributions may be useful for separating the QCD contribution from the char- 
monium resonance signal. At large invariant masses, associated production of J/?/' + 7 
should serve as a good test of hard-scattering initiated by gluons. 
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Figure captions 

Fig. 1 The scaling distribution Mf^^da/dMinv as a function of (Mj„t, — MJ)/^/s, for pp 
collisions, with beam energy = 300 GeV (solid curve) and 600 GeV (dashed curve). 

Fig. 2 Comparison of the mass distributions with data from ISR (upper figure) and 
E705 (lower figure). The theoretical curves are normalised to the data sets at 
Minv = 4.8 GeV in the upper figure, and at AM = 0.88 GeV in the lower figure. 

Fig. 3 The mass distribution da/dMinv at y/s = 23.72 GeV, as a function of Mj„i,, for pp 
collisions (left figure), and np collisions (right figure). The solid, dashed and dash- 
dotted curves correspond to cuts on the laboratory energy of the photon {E^^ > 
0, 1.0, 2.5 GeV). 

Fig. 4 The angular distribution da/dcos9*dMinv as a function of cos6'*. The curves (from 
top to bottom) correspond to Mj„i, = 3.5, 4.5, 5.5, 7.5 and 9.5 GeV. 
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